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ABSTRACT 

We use semi-analytic models and cosmological merger trees to provide the initial 
conditions for multi-merger numerical hydrodynamic simulations, and exploit these 
simulations to explore the effect of galaxy interaction and merging on star formation 
(SF). We compute numerical realisations of twelve merger trees from z = 1.5 to z = 0. 
We include the effects of the large hot gaseous halo around all galaxies, following recent 
obervations and predictions of galaxy formation models. We find that including the hot 
gaseous halo has a number of important effects. Firstly, as expected, the star formation 
rate on long timescales is increased due to cooling of the hot halo and refuelling 
of the cold gas reservoir. Secondly, we find that interactions do not always increase 
the SF in the long term. This is partially due to the orbiting galaxies transferring 
gravitational energy to the hot gaseous haloes and raising their temperature. Finally 
we hnd that the relative size of the starburst, when including the hot halo, is much 
smaller than previous studies showed. Our simulations also show that the order and 
timing of interactions are important for the evolution of a galaxy. When multiple 
galaxies interact at the same time, the SF enhancement is less than when galaxies 
interact in series. All these effects show the importance of including hot gas and 
cosmologically motivated merger trees in galaxy evolution models. 

Key words: galaxies: active, evolution, interactions, starburst, star formation - 
methods: numerical 


1 INTRODUCTION 


In the now well established model of Lambda Cold Dark 
Matter (ACDM) mergers and interactions between galax¬ 
ies play an important role in shaping their evolution 
fe.g IWhite fc Reeall978h . Both observations and theoreti¬ 


cal studies have shown that interactions between galaxies 


can c reate tidal bridges, streams and shells (agJZwicl^ 
1951 ; iToomre fc Toomra Il972l : iToomrel 1 19771 : ISmith et al 


20071). Minor mergers have been shown to t hicken disks 

'"'uinn et al.l 


I993I: Walker et al.l 

1991 

Moster et al.l I2OIOI. 

2012 


Villalobos & Helmi 

20081. 20091: 

Zolotov et al.l l20ld 

), while ma- 


jor mergers have been shown to destroy disk s and cre- 
ate elliptical and bulge dominated galaxies ilHernauistI 
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I992I: Ouinn et al.| 19931: Brook et al.l 12004: Cox et al. 

2006; 

Kazantzidis et ahlboOSl. 20091. but seelRobertson et al. 

2006; 

Governato et aLll2009l: Ouereieta et alJl2015lfor disks surviv- 


ing major mergers). 


Mergers and interactions are very common — over 70% 
of Milky Way-sized galaxies have experie n ced a 10:1 merger 
since redshift z ~ 1 IStewart et al.l (l2008l l: iLotz et ^ (1201 ll . 
e.g.). Furthermore, a significant amount of the mass formed 
in galaxies has been assembled by accreting smaller galaxies, 
varying from 10 % to >50 % of the stellar mass for galax¬ 
ies w i th Mj, > 10^°Mp) (e.g. Ose£et_S 2010j^ Lackper et al.l 


I2OI2I : iNavarro-Gonzalez et al.ll2013l l. iBell et al.l (l2006l l used 
pairs at the same redshift in the COMBO-17 photometric 
redshift survey to study the rate of interactions between 
galaxies. They found that in their sample ~ 5% of the galax¬ 
ies are in a close pair. They also found that 50% of the mas¬ 
sive galaxies have undergone a major merger since z = 0.8, 














































































2 W. Karman et al. 


and therefore conclude that major mergers play a signihcant 
role in the evolution of galaxies in this period. 

An important effect of galaxy interactions and 
merge rs is that they cause enhanced star formation 
(SFl dLarson fc Tinsle\^ll97^: Barton Gillespie et al.ll2003l : 


I Lin et aiL I I 2 OO 7 I : iHwang et al.ll201ll : |Pa.tton et al.ll2 01lh. ah 


though weak in the lo cal Universe ( Bergvall et ahl 120031 : 


iKnaDen & Jamed 

20091). and may trigger AGN activity 

fe.g. ISanders et al. 

119881: iTreister et al.ll2012l:lBessiere et al.l 


|2014) . There is increasing evidence that most extreme 
local star forming galaxies, the Ultra Luminous In- 
fraRe d Galaxies (ULIRGS), a re triggered by mergers 
fe.g. le nders fc Mirabell Il996l: iBorn^et"^ Il999l . l200d : 
iBushouse et al.l 20021 : Rupke fc VeilleuxI 2013l L and an anti¬ 
correlation is found between galaxy pair separation and the 
indicators for SF in the optical fe.g. equ i valent width of Ha , 
Sanders_e^_^ 19^B^rto^e^aU 20M; iNikolic et al.ll2004l : 


Ellison et al.ll2008r Patton et al.ll201ll L 


Galaxy mergers have been studied ex tensively usini 


Hernauisd 19921: 

Mihos & Hernauisd 

199GI: Isnringel et alJ 

2 OO 5 I: Gox et al.l 

20081: iHonkins et al. 

I 2 OO 9 IL In this ap- 


proach, two idealised (typically disk-dominated) galaxy 
models are created, each within their own idealised dark 
matter halo, and these are set up on an orbit which will 
lead to an eventual merger. The evolution during the course 
of the merger is then computed using numerical N-body and 
hydrodynamic techniques. These studies predict a strong 
link between galaxy interactions and enhanced star forma¬ 
tion. During interactions, the gravitational potential is dis¬ 
turbed, such that gas is able to lose angular momentum and 
flow inwards. Because of this, the gas density in the cen- 


place in the center of the galaxy dMihos & Hernauisd 

1996|; 

Soringel et al. 

I 2 OO 5 I: 

Cox et al.ll2008f). Cox et al.l (2008 

) and 

Hopkins et al. 

(200!; 

) quantihed this increased SF due to 


galaxy interaction and found that it is proportional to the 
mass ratio of the merging galaxies, as well as their cold gas 
fraction before the onset of the merger. Many semi-analytic 
models use the fitting formulae derived from these studies 
to model merger-driven starbursts (e.g . Bowei_e£_alJ_ 200^ 
Croton et al.ll2006l: iDe Lucia fc Blaizotll2007l : lMonaco et al l 


2007 I: ~merville et al.l 2' 


The majority of these previous idealised merger simu¬ 
lations include a cold gaseous disk in the galaxy model, but 
did not consider the hot gas component harboured in the 
dark matter halo. Full cosmological hydrodynamical simu¬ 
lations and semi-analytic models (SAMs) of galaxy forma¬ 
tion predict a large amount of hot gas around the galaxy, in 


Kauffmann et al.lll993l: iBower et al.ll2006l: ISomerville et al.l 

2008 

: Stinson et al.l 

2 OIOI. 12012: Brook et al.l 2014. but see 

also 

Cox et al.l I 2 OO 4 

for the generation of a hot halo dur- 


ing mergers). Recently, several X-ray studies found evi- 
dence for a hot g a s component around the Milky Way 
dGupta et al.l 12012 : iHodges-Kluck fc Br egmM 2013l~) and 
normal massive disc gala xies (e.g. lA~ders~ fc BregmanI 
l201ll : lAnderson et al.ll2()13l l with temperatures of more than 
10® K. 

This hot gas component was first in cluded in idealised 
hyd rodynamical merger simulations by iKim et ^ (l2009l l 
and iMoster et al.l ll201 ih . The latter found that the cool¬ 


ing of the hot gaseous halo refuels the cold gas in the center 
of galaxies, an d as a result, SF can be sustained for a longer 
ti me fsee alsolc hoi et a l.ll2014fl . Another interesting result 
of iMoster et aLn 2011 1 is that the efficiency of mergers in 
triggering SF is reduced in the presence of a hot halo. 

An additional drawback of the idealised binary merger 
approach is that initial conditions (such as the properties 
of the progenitor galaxies and the orbit) must be specihed 
a priori, and will not necessarily represent a cosmologically 
motivated distribution. In addition, such studies have been 
almost entirely restricted to binary mergers, while multiple 
mergers (a sequence of more than two mergers, in which 
the remnant does not have time to fully relax before a sub¬ 
sequent merg er) are a common ph e nomenon in the CDM 
universe (e.g. IVaisanen et al.l I 2 OO 8 I : iDuplancic et al.1 l2013l . 


MMS14). 

In t his paper we take ad vantage of the method devel¬ 
oped by IMoster et al.l (l2014h . which combines cosmologi¬ 
cal merger trees and semi-analytic models to specify well- 
motivated, cosmologically representative initial conditions 
for a set of galaxy mergers, and then carries out the evo¬ 
lution of the galaxies in the merger tree using detailed nu¬ 
merical hydrodynamics. This approach is complementary to 
fully cosmological hydrodynamic “zoom-in” simulations, in 
that very high resolution can be achieved, and controlled 
experiments to study specific physical processes can be car¬ 
ried out. This is particularly relevant for testing and rehning 
re cipes for use in semi -analytic models. We extend the work 
of iMoster et al.l ll20l4) by including the effect of supermas- 
sive black holes (BHs) and their associated feedback in our 
merger simulations and studying the modifications to the 
star formation rate (SFR) during mergers with respect to 
models without BHs. 

The paper is organised as follows: in Sec. [5] we de¬ 
scribe the methods that we employed in this work, including 
the code that has been used, the physical recipes that we 
adopted and the models that were used to create the initial 
conditions. In Sec. [3] to |4] we present our results. We first 
compare simulations with and without the hot halo compo¬ 
nent in Sec. O then we inspect the effect of multiple merg¬ 
ers in Sec. |4j and finally we present the effects of BHs in 
the simulations in Sec. [S] We discuss the performance of the 
code and possible shortcomings of the models in Sec. [S] and 
then present our conclusions from this study and possible 
improvements for future studies in Sec. 0 


2 METHODS 

In this section we briefly describe the modelling approach 
used in this paper, which is based on the simulated merger 
tree method described in detail in MMS14. The method con¬ 
sists of three main steps: 1) merger trees are extracted from 
a large dissipationless cosmological A-body simulation 2) a 
semi-analytic model is run in these merger trees to “popu¬ 
late” the halos with galaxies and predict their properties 3) 
each time two (or more) dark matter halos merge, galaxy- 
halo models are created with properties predicted by the 
SAM, placed on orbits specified by the cosmological simula¬ 
tion, and evolved forward using a numerical hydrodynamic 
code. More details about each of the three aspects are given 
below. 
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SAM 


Simulation 



dark matter 
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Figure 1. Schematic view of the scheme to combine SAMs and 
merger simulations (Fig. 6 from MMS14). On the left side the 
semi-analytic merger tree is shown with time running from top 
to bottom. At the starting time ti, a particle realisation of the 
central galaxy is created using the properties of the central galaxy 
in the SAM at this time. This system is then simulated with the 
hydrodynamical code until the time of the first merger , where a 
particle realisation of the first satellite Si galaxy is created, again 
using the properties predicted by the SAM, and introduced into 
the simulation with orbit properties taken from the cosmological 
Af-body simulation. The resulting merger is simulated until the 
next galaxies (52 and Sz) enter the main halo at t 2 , at which 
point they are also included in the simulation. This procedure is 
repeated for all mergers until the final time of the run t f. 


2.1 Cosmological N-body simulation 

The initial W-body simulation from which we extracted 
the merger trees was performed with the GADGEt2 code 
(Spri ngel 2005) assumin g WMAP3 cosmological parame¬ 
ters dSoergel et al.l [ioo^ namely; flm = 0.26, Ha = 0.74, 
h = lTo/(100 km s“^ Mpc“^) = 0.72, ag = 0.77 and 
n — 0.95. The simulation was done in a periodic box with a 
side length of 100 Mpc, and contains 512® particles with a 
particle mass of 2.8 x 10® Mq and a comoving force softening 
of 3.5 kpc. Dark matter haloes were identified in the simula¬ 
tion snapshots using a Friends of Friends (FOF) halo finder 
with a linking parameter of 6 = 0.2. Substructures inside the 
FOF groups are the n identified using the SUBFIND code 
JSpringel et aklliooil l. 


2.2 Hydrodynamical simulations 

The hydrodynamical simulations in this work were 
performed using the p arallel TREE-S PH code pgad- 
GEt3 last described by ISpringell ll2005ll . Th is code uses 
Smoothed Particle Hydrodynamics (SPH; Lucvl 1 19771 : 
lOingold Sz MonaghanI Il977l : iMonagh^ Il992 l to evolve 
the gaseous medium using an entropy conserving scheme 
llSpringel &: Hernauistl[20^ ~l. 

Important parameters in SPH simulations are the soft¬ 
ening length and the number of particles used. For the soft¬ 
ening length e we follow MMS14 and for every particle s 
species we set it to e = ei rupart /10'^° Mq llDehnenll200lll . 
where ei = 32 kpc for a particle of 10®® Mq. For one of our 


typical simulations, the number of particles are ~ 800000 
DM particles, ~ 600000 gas particles, and ~ 60000 star par¬ 
ticles, with masses of ~ 4 x 10®Mq, ~ 3 x 10®Mq, and 
~ 2 X 10® Mq respectively. In Table [AT] we give the number 
of particles for the major components in our merger trees. 

The multiphase interstellar me dium is modelled fol¬ 
lowing the SF recipe introduced by ISoringel fc HernauistI 
ll2003h . This subgrid recipe stochastically forms stars when 
cold gas clouds surpass a density threshold (pth). The newly 
formed stars eject a fraction of their material, modelling Su- 

J ierno va (SN) driven winds, following ISoringel fc HernauistI 
20031) . The star formation is thus directly coupled to SN 
feedback on a subgrid level, such that the wind mass is di¬ 
rectly proportional to the SFR, i.e. = rjMt, where rj is 
the mass loading factor which quantifies the wind efficiency. 
These gas particles receive a fixed fraction of the SN energy 
and are decoupled from hydrodynamical interactions for a 
short period of time, resulting in a constant initial wind 
speed v-w- 

The parameter s for the multiphase ISM mod el are com¬ 
puted as outlined in lSoringel fc He rnquisfl (|2003l) in order to 
match the Kennicutt-Schmidt Law I Kennicutdl 1 998l L When 
a Kroupa IMF is adopted, the mass fraction of massive stars 
is P = 0.16 resulting in a cloud evaporation parameter of 
Ao — 1250 and a SN “temperature” of Tsn = 1.25 x 10® K. 
The SF timescale is set to = 3.5 Gyr and for the galactic 
winds we adopt a mass-loading factor of t; = 1 and a wind 
speed of Vw ~ 500kms“®, which is an intermediate strength 
for winds. 

Every simulation is repeated with BHs included and 
thermal AGN feedb ack implemented as described by 
ISpringel et al.l (l2005l L BH particles ac crete mass by Ed¬ 
dington limited Bondi-Hoyle accretion (iHovle fc Lvttlet^ 
I 1939 I : iBondi fc Hovlel[l94^lBon^ll952h . and this accretion 
is assumed to deposit thermal energy into the surround¬ 
ing medium at a rate proportional to the accretion rate: 
Efeed = ££ Cr Mbh c® where Ef i s the coupling str e ngth of the 
feedback, and is adopted from iDi Matteo et akl ll2005ll . and 
Cr is the standard conversion factor of rest-mass to energy 
for a s tandard thin accretion disk from IShakura fc SunvaevI 
1119731) . 

When galaxies merge, their BH are both brought 
to the center of the remnant. During this migration, 
AGN offset from the center and pairs of AGN are pre¬ 
dicted in the remnant. Observations have shown that off¬ 
set and dual A GN are present i n a f raction of merg¬ 
ing galaxies (e.g. IComerford et akl l2009l : IKoss et I l20li : 
IComerford fc Greenel 2014h . Because the movement and or¬ 
bit of a second BH in galaxy mergers can be important for 
the evolution of the remnant, we do not reposition the BH 
particles on the local minimum of the potential. On the other 
hand we need to have stable BH particles, whose position 
is not affected by the finite resolution of the simulation (for 
example due to two body encounters with massive star par¬ 
ticles) . 

In order to overcome this problem we defined a new 
parameter named SeedRatio (SR), which defines the initial 
ratio between the seed BH mass (i.e. the actual BH mass) 
and its dynamical mass. The dynamical mass is only used by 
the gravitational computations; by increasing it with respect 
to the actual BH mass we ensure that two body encounters 
between star (or gas) particles and the BH will not sub- 
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Table 1. Summary of the parameters used for the simulations of merger trees and their fiducial values. 


Parameter 

Description 

Fiducial value 

Zi 

Redshift at the start of the simulation 

1.5 

/4-min 

Minimum dark-matter mass ratio 

0.0667/0.02 

C 

Ratio of scaleheight and scalelength of the stellar disc 

0.15 

X 

Ratio of scalelengths between gaseous and stellar disc 

1.5 

C 

Ratio of gaseous halo core radius and dark matter halo scale radius 

0.22 

dhg 

Slope parameter of gaseous halo 

0.67 

a 

Ratio of specific angular momentum between gaseous and dark halo 

4.0 

A* 

Expected final number of stellar particles in the central galaxy 

200 000 

K, 

Ratio of dark matter and stellar particle mass 

15.0 

-^reSjSat 

Ratio of satellite and central galaxy particle mass 

1.0 

Amin 

Minimum number of particles in one component 

100 

ei 

Softening length in kpc for particle of mass m = 10^® Mq 

32.0 

f* 

Gas consumption time-scale in Gyr for SF model 

3.5t 

Ao 

Gloud evaporation parameter for SF model 

1250.0t 

dsF 

Mass fraction of massive stars for SF model 

0.16t 

Tsn 

Effective supernova temperature in K for feedback model 

1.25 X 10*t 

V 

Mass loading factor for wind model 

1.0 

"^^wind 

Initial wind velocity in kms~^ for wind model 

500.0 

SeedRatio 

Ratio of seed mass and dynamical mass of BH particle 

0.01 


t The SF parameters assume a Kroupa IMF. 

stantially modify the BH trajectory. On the other hand all 
physical processes associated with the BH (e.g. feedback) 
are governed by its real mass and are not influenced by our 
choice of the initial SR ratio. 

It should be noted that as the BH gradually accretes 
mass at a subgrid level, only the actual mass of the BH 
increases, while when another particle is swallowed by the 
BH particle the dynamical mass increases. We found that a 
ratio of 1000 for dynamical masses of BH particles to gas 
particles ensures very stable BHs orbits in all our merger 
simulations. At our resolution this results in SR = 0.01 
(meaning the BH dynamical mass is 100 times the BH mass). 
An overview of parameter values used in this work is given 
in Table [T] 

Before moving to the next section, we would like to com¬ 
ment on the possible effects of the shortcomings of the pgad- 
GEt3 im plementation of the S PH method on our results. 
Recently, iHavward et al.l ll20l4l showed that the agreement 
between the adaptive mesh refinement (AMR) approach and 
SPH methods is very good with regard to SFRs in galaxy 
mergers. They also show that the inclusion of BHs does not 
significantly modify the small difference in the SFH in SPH 
and AMR methods. This shows that the conclusions in this 
paper should be robust to the numerical treatment of hy¬ 
drodynamics, in spite of recent concerns about the accuracy 
of the classical SP H formulation in certain situations (e.g 
lAgertz et al.ll2nn^ . 

2.3 Semi-analytic models and Simulated Merger 
trees 

In thi s paper we use the same approach as iMoster et al.l 
(120111) and MMS14, and refer to those papers for a more ex¬ 
tensive discussion. In our method, presented schematically 
in Fig. [T] we build a dark matter tree fr om a large-scale 
Ai-body simulation and use the SAM of ISomerville et ah] 


ll2008l) to populate the A-body merger tree with galaxies. 
The SAM predicts the baryonic properties for each galaxy 
within the halo merger tree by applying simple but physi¬ 
cally motivated recipes for cooling and accretion, star for¬ 
mation and chemical enrichment, stellar feedback, black 
hole growth and AGN feedback, and morphological trans¬ 
formation and starbursts in mergers. The radial sizes of 
disks are modelled using an angular momentum based ap¬ 
proach which has been shown to predict a size-mass rela¬ 
tion th at is in very good agree ment with observations up to 
z ^ 2 dSomerville et all 120081 ). These SAMs have been ex¬ 
tensively tested and shown to reproduce many fundamental 
statistical properties of galaxies, including stellar mass func¬ 
tions and luminosity functions from z ~ 0-6, disk gas frac¬ 
tions, and the fraction of different morphological types from 
z ^ 0-3 dSornerville et all l2008l . l20ld : IPorter et al.l l20l4 
iBrennan et al.ll2015l ). 

In Fig. m we show an example of such a merger tree, 
where time runs from top to bottom. At a chosen time ti we 
extract a halo with its baryonic properties and merger tree. 
In this simple example the main system experiences four 
mergers, at H, t 2 (two) and At ti we obtain the proper¬ 
ties from the SAM, and create a particle realisation of the 
main halo (Ci), represented by the brown arrow. In Fig. [T] 
this generated galaxy is shown in the top right, and consists 
of a dark matter halo (grey), a hot gaseous halo (red), a cold 
gaseous disc (blue), a stellar disc (yellow) and a stellar bulge 
(green). We then simulate the evolution of this galaxy up to 
ti, at which point we add a particle realisation of satellite 
galaxy Si to the simulations. While we obtain all baryonic 
properties from the SAM, the orbital parameters of Si are 
directly taken from the dark matter A-body simulation. In 
this way we create a cosmologically motivated merger his¬ 
tory. 

After the inclusion of this satellite galaxy, we evolve the 
system in our simuations up to t 2 . At t 2 we add two more 
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satellites {S 2 and S3) at the virial radius using the same 
method as before and continue the evolution of the system. 
We repeat this procedure up to the end of the simulation 
at tf. By using this method we naturally include multiple 
mergers, such as at t 2 and at ts where one of the satellites is 
still orbiting when satellite S4 enters the halo. The growth of 
the main halo occurs through two channels. First, the main 
halo accretes gas and dark matter, which takes place from ti 
to ti or from ts to tf in the given example. Second, smaller 
satellite galaxies feed the main halo through mergers, e.g. 
from ti to t 2 . 


2.3.1 Galaxy Initial conditions 

We follow the descrip tion provided by Springel et al.l ll2005li , 
with the extension of lMoster et al.l (l201ll ) to create particle 
realisations for our initial conditions. Each galaxy consists 
of two stellar components (a disc and a bulge with masses 
Mdiscand Mb), two gaseous components (a cold disc and a 
hot halo with masses Meg and Mhg), and a dark matter halo 
with mass Mdm. 

The stellar and gaseous discs are set up such that they 
are rotationally supported and have surface density profiles 
that decline exponentially outwards. The parameter x cou¬ 
ples the scale length of the cold gaseous disc to that of the 
stellar disc by Vg = xrd, and the vertical structure of the 
discs is determined by a radially independent sech^ profile 
with scale height 20 . The velocity dispersion in the stellar 
disc is set equal in the vertical and radial direction, while 
the gas temperature is determined by the equation of state 
(EOS). A balance of the pressure due to the EOS and the 
gravitational potential then determines the vertical struc¬ 
ture of the disc selfconsistently. 

For the stellar bulge and dark matter halo, we assume 
a lHernauisd lll990l i profile with scale lengths rs and Ts. The 
dark matter halo is characterised by a concentration param¬ 
eter c = rvir/^s and a halo spin A. We assume that the stellar 
bul ge is spherical and n on-rotating. 

iMoster et al.l (l201ll ') calibrated the properties for a hot 
gaseous halo in order to match the cosmic star formation 
rate. We adopt a spherical slowly-rotating density profile, 
that follows the observationally motivated /3-pr ofile (e.g. 
ICavaliere fc Fusco-Femianol Il97d : lEke et al]Il998li . Requir¬ 
ing that pressure supports an isotropic model and hydro¬ 
static equilibrium within the gravitational potential fixes the 
temperature profile of the hot gaseous halo. A further cal¬ 
ibration needed was the spe cific angular momen tum of the 
gaseous halo, and we follow iMoster et al.l (l201ll l in setting 
the specific angular momentum of the gaseous halo ji'hg to 
a multiple a of the dark-matter specific angular momentum 
jdm, such that jhg = ajdm, with a = i. 

In simulations with BHs, the initial mass of the BHs is 
extracted from the SAM, but due to resolution limitations 
we only include BHs with a mass larger than 10®M©. All 
the stellar and gaseous parameters (e.g. Mdiac, Mb, Meg, 
Mhg, Vd, etc.) are taken from the SAM predictions, while 
the dark matter parameters (Mdm, c, A) and all the orbital 
parameters (position and velocity) of the mergers are taken 
directly from the A-body simulation. 


2.4 Merger tree selection 

We set the lower mass limit of satellite galaxies to be in¬ 
cluded in the numerical realisation to 1/15 of the mass of 
the central galaxy. We chose this limit because it has been 
previously shown that mergers with a mass ratio > 5 

barely influence the SFR lICox et al.l 1200811^ Although we 
could therefore exclude any satellites with > 5 from 

our simulations, we choose to be conservative and set the 
limit at ~ 15, in order to catch any possible secondary 

effects on the merger trees of slightly smaller satellites. 

From the A-body-fSAM simulation we selected 12 
merger trees based on the following criteria: 

• The final mass of the halo has to be of the same order 
of magnitude as that of the Milky Way (e.g. 2 — 3 x IO^^Mq) 

• The merger tree must contain at least one major merger 

• The total sample has to contain mergers of every mass 
ratio from 1 to 5 in steps of 0.5. 

• The galaxy with the smallest mass ratio should enter 
the virial radius after redshift 1.5 


Milky Way mass galaxies sit at the peak of SF 
efficiency (Fig. [^, and the transition between disc to 
spheroid domina ted systems also occurs at this mass scale 
llConselicel l2006f l. Both of these effects are strongly re¬ 
lated with mergers. In addition, the steep relation between 
halo and stellar mass at lower masses results in a drop 
of the influence of major stellar mass mergers for lower 
mass galaxies. At higher mass, it has been shown that 
additional f eedback recipes are needed to prevent over¬ 
cooling (e.g. Venmleo&ke^nol^ 20061: Gabor et al.l I 2 OIII : 


IVoeelsberger et ^ 20141: Schave et al.l 201^ 


We invoke the second and fourth criteria to satisfy the 
aims of our study, i.e. determine the effect of the merger ratio 
on the star formation history in cosmologically motivated 
mergers at 2 < 1.5. The third criterion ensures that we 
sample the full range of merger ratios. 

The final merging process of the two most massive 
galaxies is not disturbed significantly by other galaxies in 
the majority of the trees, so during the final merger they 
are similar to binary mergers. However, we selected three 
merger trees that harbor multiple simultaneous mergers, in 
orde r to study th e importance of three body interactions 
(e.g. lMorendl2012l '). 

We also evolved the central galaxy in isolation for the 
same time interval, such that we have a control sample with¬ 
out mergers. We will refer to these simulations as the “iso¬ 
lated” cases. Each merger tree is simulated from a starting 
redshift 2 = 1.5 down to 2 = 0. Finally we reran all merger 
trees with the addition of BHs and AGN feedback. All to¬ 
gether we have a suite of nearly 100 different galaxies and 
simulations, leading to a substantial statistical sample. 


3 STAR FORMATION EVOLUTION 

We computed the stellar-to-halo-mass ratio for all of our 
initial and final central galaxies, and compared it to the 


^ We did run three merger trees with a cut of 1/50 to validate 
this assumption and found no substantial differences in the SF 
history. 
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Tree nr 

f^host 

W^°Mq 

Msat 

10^° Mq 

h 

^*,host 

^*,sat 

661wsl 

12.7 

2.04 

6.23 

0.385 

0.0895 

661ws2 

24.7 

8.39 

2.95 

2.18 

0.955 

661bhl 

13.9 

2.11 

6.59 

0.313 

0.226 

661bh2 

22.7 

7.85 

2.89 

1.81 

1.05 

872ws 

27.4 

9.01 

3.04 

4.72 

2.95 

872bh 

23.2 

6.47 

3.59 

3.86 

2.01 

990ws 

21.6 

8.98 

2.41 

2.85 

2.10 

990bh 

19.6 

8.80 

2.22 

2.42 

2.13 

1166wsla 

3.21 

0.77 

4.17 

0.781 

0.104 

1166wslb 

3.21 

1.93 

1.66 

0.781 

0.577 

1166bhla 

3.46 

1.05 

3.30 

0.650 

0.132 

1166bhlb 

3.46 

2.13 

1.62 

0.650 

0.494 

1178WS 

7.42 

2.66 

2.79 

0.700 

0.452 

1178bh 

7.46 

2.88 

2.59 

0.657 

0.320 

1188wsl 

5.14 

3.77 

1.36 

0.914 

0.508 

1188ws2 

9.68 

1.27 

7.62 

2.83 

0.365 

1188ws3 

11.7 

7.41 

1.58 

3.63 

1.66 

1188bhl 

4.85 

3.49 

1.39 

0.825 

0.487 

1188bh2 

7.72 

0.948 

8.14 

2.25 

0.364 

1188bh3 

9.15 

4.07 

2.25 

2.77 

6.60 

1415wsla 

1.04 

1.19 

0.87 

0.197 

0.221 

1415wslb 

1.04 

1.50 

0.69 

0.197 

0.269 

1415ws2 

7.13 

4.71 

1.51 

1.95 

1.05 

1684ws 

10.8 

5.33 

2.02 

2.78 

0.605 

1684bh 

11.1 

4.79 

2.32 

2.81 

0.567 

2096ws 

4.49 

3.95 

1.14 

0.728 

0.373 

2096bh 

4.36 

3.34 

1.31 

0.675 

0.352 

3747wsl 

4.11 

0.605 

6.79 

0.182 

0.0414 

3747ws2 

3.97 

2.82 

1.41 

1.39 

0.366 

3747bhl 

3.72 

0.856 

4.35 

0.167 

0.106 

3747bh2 

3.83 

2.34 

1.64 

1.28 

0.199 


Table 2. Properties of the progenitor galaxies immediately before 
they merge. Col. (1): the id number of the merger tree, (2) total 
mass of host galaxy within iJghosti (3) total mass of satellite 
galaxy within i?s,sat, (4) merger ratio /r = , (5) stellar mass 

of host within i?s,host i (6) stellar mass of satellite within i?s,sat. In 
Col. 1, the label ws means that the simulation was run with only 
the hot halo model, while the bh label means a simulation with 
hot halo and BHs. Numbers after these labels give the number of 
the merging event, and letters a subset of multiple merger events. 
Because merger trees 780 and 2809 do not have any mergers by 
^ = 0, they are not included in this table. 


constraints derived bv iMoster et al.l ll2013l l using abundance 
matching. This comparison shows that the initial conditions 
and final properties of our galaxies are consistent with these 
constraints. 

The results are shown in Fig. [2] although the initial 
conditions are mostly located within la of the relation, the 
galaxies at redshift z — 0 mostly lie above the Ict region. 
This is not unexpected because we selected merger trees 
with active merger histories, which will have preferentially 
larger amounts of mass growth and SF than the average. 
Merger trees run with BHs do agree better with the stellar- 
halo mass relation, as feedback from AGNs reduces the SF. 
Overall our galaxies do not seem to suffer from a serious 
overcooling problem and are consistent with observations. 

We simulated a total of twelve merger trees from red- 
shift z = 1.5 to 2 = 0, and summarise the characteristics of 
all mergers in Table [2] The merger ratio ji = is an 

important parameter in this study. Satellite galaxies, after 
entering the halo of the central galaxy and before merging 


with it, lose m ass due to dynamical friction and ram pressure 
stripping (e.g. Bovlan-Kolchin et al.ll2008l : IWeinmann et all 
l20ia IChang et al.l l2013ll . This mass evolution introduces 
some arbitrariness in the definition of the merger ratio, be¬ 
cause it depends on when it is computed. We decided to 
compute it right before the merger between the galaxies, in 
order to take into account the mass loss. We first fitted the 
dark matter pro file of both galaxies, before the merger, to 
an NFW profile (iNavarro et al.lll997h and obtained a value 
for the scale radius r^. We then calculated the total mass 
(dark matter-|-gas-|-stars) inside Vs at the time of merging 
(ti) for both satellite and central, and used these masses to 
define p. 

The mass ratios obtained with this method are in very 
good agr eement with results from substructure finders like 
SUBFIND dSpringel et al.l[2001^ . This method was always able 
to distinguish between the central and satellite galaxies at 
very close encounters, whereas this proved difficult for SUB- 
FIND. 


3.1 The effect of the Hot Halo 

As a first test, we compare simulated merger trees with and 
without the inclusion of the (observationally expected) hot 
gaseous halo component. We present the SFR for merger 
tree 990 in Fig|3] for a run without hot halo (left panel) and 
with hot halo (right panel). In each panel the results of the 
merger simulation are shown in red, while the isolated run 
(e.g. central galaxy without satellites) is shown in green. It 
is important to note that the merger ratio of two galaxies 
can vary after the removal of the hot halo, due to a different 
fraction of hot gas. 

The cold-gas- only simulation i s in agreement with pre¬ 
vious studies (e.g. ICox et ^l2008li : the SFR drops rapidly 
when the gas is depleted, and a starburst occurs during the 
final merging process. The starburst increases the SFR by 
an order of magnitude and there is almost no SF after the 
merger, because all the gas has been consumed during the 
starburst. The quenching of SF due to starvation is partic¬ 
ularly clear in the isolated run, where we see a progressive 
fading of SF as time goes by. We caution that the relative 
increase in S F can be less than an order of a magnitude, 
for example, iDi Matteo et al.) (l2007l . l2008fl find an average 
increase of a factor of ~ 3 during mergers. 

The simulation including a hot halo (right panel) shows 
a qui te different behavior (as already noted bv iMoster et al.l 
I 2 OIII) . The hot halo cools and therefore the cold gas in the 
center of the galaxy is refuelled and SF can be sustained 
for the whole duration fo the simulation; this is particularly 
clear in the isolated run. 

Similar to MMS14, we see an important difference be¬ 
tween the simulations with and without a hot halo during 
the merger. The starburst occurring in the simulation with a 
hot halo produces many more stars than the one in the sim¬ 
ulations without a hot halo, however, the relative increase 
with respect to the isolated galaxies is smaller. In the simu¬ 
lations without a hot halo we observe an order of magnitude 
difference between the isolated and merging galaxies, which 
is significantly larger than the factor of ~ 6 difference in the 
simulations with a hot halo. This difference is mostly due to 
the different evolution of the SFR in the isolated runs. The 
isolated galaxy without the hot halo depletes a large frac- 
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Halo mass [M/Mg] Halo mass [M/Mg] 

Figure 2. Stellar mass versus dark-matter mass within the virial radius. The left plot shows the initial galaxy properties obtained from 
the SAM, while the right plot shows the resulting simulated central galaxies at redshift zero. The red stars show the simulations with the 
hot halo model and t he blue circles show the simulations with the hot halo model and BHs. Overplotted are the abundance matching 
constraints derived bv iMoster et ahl ll201.'ll ~). The dark shaded area shows the la range on the constraints, and the light shaded area is the 
2(7 range. 


tion of its gas by forming stars in the first 2 Gyrs (before the 
merger happens), but when it has a hot halo surrounding it, 
it is able to keep forming stars for a longer period (because it 
replenishes its cold gas reservoir). This reduces the relative 
efficiency of a merger in triggering additional SF. 

In Fig. 0] the SF History (SFH) of merger tree 780 is 
shown with and without the hot halo component. For the 
simulation without a hot halo, the behavior is very simi¬ 
lar to tree 990, with an increased SF at t ~ 10 Gyr. This 
“bump” in the SF clearly shows the effect of the merger, 
but it is not present in the simulation with a hot halo even 
though the mass ratio is very similar. When looking directly 
at the simulation snapshots at ^ = 0, we see indeed that the 
two galaxies did not merge and are at a distance of 73 kpc, 
despite two close passages (with distances of -^100 kpc and 
<5 kpc) and in contrast to expectations based on the dark 
matter only simulation. 0 

Even more interesting is the decrease of SF in the merg¬ 
ing hot halo simulation compared to the isolated galaxy. 
This decrease is seen over the full simulated period. This 
result is in contrast to the findings of previous m erger simu¬ 
lations (e.g. ICox et ^l2008l : lHopkins et al.ll2009l 'l. where the 
SFR only experienced a positive boost due to interactions 
of galaxies. In order to better understand the origin of this 
SFR reduction we looked in more detail into the cold/hot 
gas fraction in the merger simula t ion. _ 

ISinha fc Hollev-BockelmannI (l2009f ] and iMoster et al.l 
||201iF already found that when two galaxies merge the hot 
halo component experiences heating due to the creation of 
shocks and the dissipation of the satellite galaxy orbital en¬ 
ergy. This results in a net increase of the gas temperature 
and in a longer cooling time. 

In Fig. [S] (bottom right panel) we show the tempera¬ 
ture for the halo gas at a distance from 30 to 40 kpc from 
the galaxy center. This temperature is shown for the central 


galaxy in isolation (green) during the merger (red) and for 
the incoming satellite galaxy again in isolation (blue) and 
in the merger simulation (cyan). For both galaxies (central 
and satellites) the gas temperature is higher i n the merger 
run th an in isolation, confirming the results of iMoster et al.l 
ll201lll . This higher gas temperature increases the cooling 
time, and reduces the amount of cold gas available for SF, 
as shown in the bottom left panel of Fig. [S] For the satel¬ 
lite galaxies the effects of ram pressure and tidal stripping 
make the difference between the isolated and merger cases 
even more pronounced. Finally the increased temperature of 
the hot gas and the subsequent reduced amount of cold gas 
decreases the SFR and produces a lower final stellar mass 
(right and left upper panels in Fig. in the merger run with 
respect to the isolated case, making the merger an overall 
“inhibitor” for SF. 

The increase of the SFR in the isolated host galaxy at 
~ 6 Gyr is caused by disk instabilities that originate from 
the accretion of cooling gas. Because the cooling rate is lower 
in the merger, the amount of accreted gas is insufficient to 
introduce the instabilities observed in the isolated case after 

6GyiE|. 

Another important effect that might be at play in this 
galaxy is tidal stripping. If the first encounter at 4.8 Gyr 
generates strong tidal effects, a significant fraction of the 
cold gas can be expelled from the galaxies. This stripping 
would make the gas unavailable for star formation, and pro¬ 
vides another method of reducing star formation that pro¬ 
posed above. However, when visually inspecting this galaxy, 
we find no signs of strong tidal effects, such as tails, during 
the first encounter. This shows that the difference in star 
formation between the isolated and merging galaxies is not 
dominated by tidal effects. 

When looking at the complete suite of 12 merger trees, 
we see a number of important trends. First, in all simula- 


^ Since in this merger tree, and also in merger tree 2809, the two 
galaxies do not merge by z = 0, we did not use these merger trees 
in the remainder of the paper. 


® The distance between the two galaxies increases from 200 kpc 
at 5.5 Gyr to 500 kpc at 7.5 Gyr, showing that during this period 
there is no significant interaction. 
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Figure 3. The SFR of merger tree 990. The left panel shows the simulations where no hot halo is present, and the right panel shows the 
simulations with a hot gaseous halo around the galaxy. The red line shows the SFR for the primary galaxy with the simulated merger, 
and the green line is the primary galaxy simulated in isolation. The vertical dashed orange line indicates the time when the merger occurs 
while the blue dotted line shows flyby’s. 




Figure 4. Same as Fig. [3]but for merger tree number 780. 


tions the inclusion of the hot halo component increases the 
overall amount of SF both in the isolated runs and during 
mergers. The hot halo is able to refuel the cold gas reser¬ 
voir and prevent the artificial quenching of SF due to gas 
depletion. Second, starbursts of varying sizes always occur 
during close encounters between galaxies. Mergers are able 
to funnel gas into the central region of the galaxy, where it is 
rapidly converted into stars. Most starbursts show a factor 
of ~ 6 increase in SFR with respect to the isolated run of 
the same galaxy (see next section for a more quantitative re¬ 
sults). Finally we found that mergers can also decrease the 
total SFR especially over long periods. Orbiting satellites 
increase the temperature of hot halo gas, by both the cre¬ 
ation of shocks and transfer of orbital energy into thermal 
and kinetic energy of the gas through dynamical friction and 
hydrodynamic drag. The increased temperature of the hot 


halo leads to longer cooling times, diminishing the amount 
of hot gas able to cool and reducing the overall amount of 
SF. 


3.2 Starburst efficiency with hot halo 

In MMS14 it was shown that star formation enhancement 
in multiple mergers can be significantly different than in a 
series of binary mergers. Here we extend the work of MMS14 
on this issue. In this work we quantify our results in a 
more useful and descriptive way, so as to better compare 
the star formation during mergers in isolated binary merg¬ 
ers and in multiple mergers as predicte d by cosmo l ogica l 
mer ger trees. Foll o wing the approach of I Cox et ^ (l2008fl 
and iMoster et all (1201 il l . we introduce the starburst effi¬ 
ciency parameter esurst, dehned as the fraction of cold gas 
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Figure 5. Merger tree 780: the top left plot shows the evolution of the SFR of this merger tree from redshift 1.5 to 0; the top right panel 
is the stellar mass over this period; the bottom left shows the amount of cold gas present in the galaxy, and the bottom right panel is the 
average temperature of the gas between 30 and 40 kpc around the galaxy. The red (green) line always shows the primary (secondary) 
galaxy in the merger simulation, while the blue (cyan) lines shows quantities for the primary (secondary) galaxy simulated in isolation. 


present in the pre-merger galaxy that is consumed by the 
additional burst of SF. 


_ Afburst 

Cburst — -TT 

premerger, cold, gas 


( 1 ) 


Where the burst mass Mburst is defined as the additional 
amount of stellar mass formed due to the merger between 
and t 2 , which defines the beginning and the end of the burst 
(e.g. ti = 7.5 and t 2 = 8.5 in the left panel of Fig. [3l and 
should not be confused with and t 2 from Fig.[T] which are 
the entry times of different satellites into the simulation): 


Mburst = ^ SFi?(t)dt-i [SFR{ti) + SFR{t2)] (tz-ti), (2) 


The burst masses and burst efficiencies from our suite of 
mergers are listed in Table [3] In addition, this table shows 
the bulge fraction /bulge, the gas fraction /gas, and other 
important properties of our mergers. The bulge fraction is 
calculated as described in detail by Kannan et al. (in prep.), 
but is based on a decomposition of the stellar component in 
a bulge and a disk, using the angular moment of the stel¬ 
lar pa rticles, similarly to the approach of lScannapieco et al.l 
ll201llf . Fig.[6]shows the relation between the burst efficiency 
and the galaxy merger ratio /i. For /r < 5 there is a large 


scatter in the values of eburst, while for p > 5 the starburst 
efficiency seems to converge to values around 0.15 (though 
the scatter is difficult to assess due to small number statis¬ 
tics). 

Following ICox et al.l (120081') . we performed a power law 
ht to the burst efficiency as a function of mass ratio: 



We find that for our simulations including the hot halo the 
best fit is found for em = 0.52 and 7 = 0.81, and the result¬ 
ing relation is shown in Fi g. [6l along with t he ht with the 
original values proposed by Cox et all ll2008l . Eq. 5). 

The valu e s for e in and 7 are comparable to those found 
bv ICox et al.l ll2008h . but the large scatter in the eburst val¬ 
ues suggests that other parameters besides the mass ratio are 
impacting the burst efficiency. The larger scatter is due to 
three factors. First, t he number o f mer gers included in this 
study is larger than in ICox et al.l (l2008l l. who ht to ten data 
points. Second, we use a range of cos mologically based orbits 
and initial galaxy prop erties, whe reas Cox et ahl (l2008f) only 
varied the mass ratio. ICox et ^ (l2008l f show that the gas 
fraction, bulge-to-disk ratio and pericentric distance also in- 
huence the starburst efficiency, however, this is insufficient 
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Tree nr 

ti 

Gyr 

t2 

Gyr 

h 

fgas 

fbulge,host 

^bulge,sat 

■l^burst 

Mcold 

IOIOM0 

^burst 

661 

5.23 

6.34 

6.23 

0.65 

0.001 

0.16 

0.232 

1.933 

0.12 

661 

7.87 

9.54 

2.95 

0.32 

0.20 

0.001 

0.135 

3.631 

0.037 

872 

9.78 

10.51 

3.04 

0.10 

0.23 

0.56 

0.171 

2.645 

0.064 

990 

7.26 

8.47 

2.41 

0.24 

0.19 

0.18 

1.079 

2.056 

0.52 

1166 

7.60 

8.30 

4.17 

0.50 

0.22 

0.33 

0.123 

1.161 

0.11 

1166 

8.30 

8.85 

1.66 

0.26 

0.43 

0.11 

0.0365 

1.444 

0.025 

1166t 

6.97 

8.52 

1.12 

0.52 

0.005 

0.008 

0.490 

1.407 

0.35 

1178 

7.80 

9.80 

2.79 

0.52 

0.047 

0.18 

0.538 

2.854 

0.19 

1188 

6.38 

7.59 

1.36 

0.44 

0.061 

0.079 

0.708 

1.923 

0.37 

1188t 

6.41 

7.36 

1.52 

0.44 

0.10 

0.10 

0.736 

1.430 

0.51 

1188 

10.27 

10.76 

7.62 

0.22 

0.24 

0.15 

0.0617 

1.141 

0.054 

1188 

12.85 

13.72 

1.58 

0.16 

0.42 

0.14 

0.172 

1.923 

0.089 

1415tt 

7.70 

8.85 

0.87 

0.55 

0.26 

0.08 

0.671 

1.139 

0.59 

1415tt 

7.70 

8.85 

0.69 

0.55 

0.26 

0.11 

0.671 

1.139 

0.59 

1415 

11.65 

13.00 

1.51 

0.18 

0.45 

0.10 

0.458 

0.804 

0.57 

1684 

8.33 

11.08 

2.02 

0.33 

0.15 

0.002 

1.012 

3.120 

0.32 

2096 

7.81 

10.03 

1.14 

0.44 

0.055 

0.024 

2.125 

2.838 

0.75 

3747 

5.25 

6.45 

6.79 

0.81 

0.009 

0.013 

0.247 

1.525 

0.16 

3747 

9.95 

11.40 

1.41 

0.37 

0.21 

0.20 

0.813 

3.058 

0.27 


Table 3. Results of the hot halo only simulations, where only ten out of the initial twelve merger trees contained major mergers. The 
first column gives the tree number, the second and third the begin and end of the starburst, the fourth the mass ratio from Table [2] In 
column five, six, and seven the cold-gas fraction within i?s of both galaxies, the bulge fraction of the host, and the bulge fraction of the 
satellite are given. The eighth column is the burst mass from equation [2] and the ninth is the amount of cold gas in the central region. 
The burst efficiency as defined in equation [T] is given in the last column, (t): The least massive of the three simultaneously interacting 
galaxies is removed and the merger is reran, see Sec.|2 (^^)' Combination of three galaxies, unable to distinguish contributions, so both 
merger ratios are given, see also Sec. 


to explain the large scatter. This is most clearly demon¬ 
strated by comparing merger trees 990 and 1178, which have 
similar properties, but significantly different burst efficien¬ 
cies. Third and possibly most important, the more complex 
setting of the full merger tree influences the state of the 
galaxies. 


In Fig. [7] we examine the origin of the scatter in more 
detail, and we show that the departure from the fitted power 
law cannot be explained by a single parameter. The scatter 
around the power law shows only weak trends with /gas, 
or /bulge- We use the duration of strong interaction, At, as 
a tracer of orbital influence, rather than trying to describe 
the orbit by an ellipse. If one did describe the orbit with 
an ellipse, the parameters of the ellipse would change be¬ 
tween every time step due to interactions with other galax¬ 
ies. Therefore, the initial orbital parameters ar e a very poor 
trace r of the final orbital parameters (see also iTsatsi et al.l 
l2015h . The final orbital parameters are better traced by At, 
where a more tangential orbit has a higher value of At. Al¬ 
though the overall distribution of the difference in burst 
efficiency seems random, there might be a trace of a pos¬ 
itive trend with increasing At when we only select mergers 
with a similar merger ratio (1 < /r < 2, large points in Fig. 
(Tjl. The positive correlation between the duration of the en- 
counter and the s treng th of the starburst was also found by 
iDi Matteo et ahl (120071 '). In summary, this shows that many 
factors contribute to the efficiency of the merger-triggered 
starburst and it may be difficult to devise a simple fitting 
formula. We will show in Sec. [5] that this conclusion holds 
after the addition of BH feedback to the simulations. 


Tree 

Gyr 

t2 

Gyr 

h 

-l^burst 

■l^cold 

^burst 

780 

9.35 

11.65 

1.03 

0.111 

0.340 

0.328 

872 

11.00 

12.60 

2.99 

0.008 

0.341 

0.023 

990 

7.55 

9.70 

2.43 

0.187 

1.066 

0.176 


Table 4. Results of the simulations without a hot halo. Only 
three out of twelve merger trees were simulated with these initial 
conditions. The columns correspond to columns 1, 2, 3, 4, 8, 9, 
and 10 from table [3] 



Figure 6. Mass ratio versus starburst efficiency. The red points 
are the mergers in the simulations with a hot halo, while the green 
squares show the simulations without a hot halo . The black line 
is the starburst efficiency relation of ICox et al.l ll2008l l and the 
green line represents the least squares fit to our results. 





















Star formation during mergers 


11 




Figure 8. The multiple mergers in tree 1166. The top left plot shows the SFR, the top right plot the total stellar mass, and the bottom 
left plot the cold-gas mass in the central volume with radius of 20 kpc. The bottom right plot shows the average temperature of the halo 
gas between 30 and 40 kpc (T 30 ). The red line is the central galaxy in the original simulation, and the green line shows the simulation 
where only the central galaxy and first satellite is included. The blue line shows the simulation where the first satellite is left out, and 
only the two most massive galaxies are included (the central galaxy and SAT2). The vertical orange dash-dotted line and pink dashed 
line indicate the time when the two satellites finally merge with the host galaxy (p = 4.17 and p = 1.66) in the simulation with a triple 
merger, while the vertical dotted blue line indicates the time when the galaxies merged in the binary merger simulation including only 
the most massive galaxies (/i. = 1.12). 
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Figure 7. Difference between the burst efficiency compared to 
the power law fit to our results shown in Fig. O Aeturst = 
^burst ~ ^burst PL- The difference is shown as a function of the 
gas fraction /gas, the bulge fraction of the main galaxy /bulge, 
and the duration of the starburst At = t 2 — ti, which is a tracer 
of the orbital properties. The large circled points are all major 
mergers with 1 < p < 2, while the smaller points are all more 
minor mergers. 


4 MULTIPLE MERGERS 

In addition to situations where only two galaxies interact si¬ 
multaneously, so-called binary mergers, it also happens that 
more than two galaxies interact with each other or share a 
common halo. These multiple merger events are com mon in 
the LCDM universe (e.g., see IVaisanen et al.ll2008l . for an 
observed merging triplet), and can have important conse¬ 
quences for galaxy evolution. In more than 50% of all merger 
pairs since « < 5, a second satellite enters the common halo 
within five dynamical times after the first satellite entrance 
(MMS14). 

Nearly all previous studies on the effect of major merg¬ 
ers on SF focused only on binary mergers. Here we use our 
cosmologically motivated merger tree to study the effects 
of multiple strongly-interacting mergers on SF. We selected 
three merger trees (namely 1166, 1188, 1415) where at least 
three galaxies strongly interact with each other during a 
merger. We refer to these multiple strong interactions as 
multiple merger events. 
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In Fig. [8] we present the results for tree 1166. In each 
panel three different simulations are shown: 1) the multiple 
merger where the central galaxy first merges with the first 
satellite (SATl) and then the second satellite (SAT2) enters 
the virial radius before the end of the first merger; 2) the 
binary merger between the central galaxy and the first satel¬ 
lite (without SAT2); 3) a binary merger between the central 
galaxy and SAT2, which happens to be more massive than 
SATl. 

The effect on SF is stronger in the merger of the two 
more massive galaxies (central and SAT2) than in the case 
of the triple merger. In the binary merger the starburst is 
more pronounced and more extended. The first satellite has 
a very minor impact on the SFR, moreover in the case of 
a triple merger SATl has an overall negative effect on the 
SFR, which is reduced with respect to the binary merger 
between central and SATl. 

Although the stellar mass formed in the starburst of 
the merger with only SAT2 is larger than in the original 
simulation, the stellar mass that resides in the galaxy at the 
end of both simulations is very similar. The difference in 
stellar mass formed is compensated by the accretion of the 
stellar mass native to SATl, and therefore results In similar 
total stellar masses. 

The reduced SF in the multiple merger is clearly due to 
the reduced amount of cold gas, as shown by the third panel 
in Fig. [8] In the multiple merger simulation the orbiting of 
the third galaxy increases the temperature of the hot gas 
(see the fourth panel of Fig.[^ delaying the cooling process. 
In Fig. I^we show the distribution of gas in this merger tree 
at important stages of this interaction. 

Fig. llUl shows the same panels as Fig.[S]for merger trees 
1188 (upper row) and 1415 (lower row). In these plots we 
see similar results. In 1188 SATl is the most massive and 
indeed produces a larger starburst than the triple merger, 
in agreement to what we found in 1166. In merger tree 1415 
we encountered a more complex situation. From Table [5] 
it can be seen that the merger ratio between the central 
galaxy and the two satellites SATl and SAT2 is less than 
1 using our definition of fi. This means that both the satel¬ 
lites are more massive than the central galaxy, and that the 
merger between the two satellites will be triggering the most 
SF, also in the triple merger. Simply replacing the central 
galaxy by one of the satellites and keeping the same orbital 
parameters, does not result in a merger before 2 = 0. The 
orbital parameters would therefore have to be significantly 
adjusted in order to have a merging binary system of the two 
satellites, and this would complicate the comparison. For all 
these reasons it is very difficult to compare the results of 
1415 with the other multiple mergers. 


5 EFFECT OF BLACK HOLE FEEDBACK 

Feedback from accretion supermassive black holes can 
hav e an important effect in regulating SF during merg- 


ers ( 

Di Matteo et al.ll2005l: iHonkins et al.ll2006l: iGhoi et al.l 

|201i 

). During mergers, gas can be efficiently funnelled to- 


wards the center of the galaxy triggering rapid accretion 
onto the BH (Quasar phase) and a subsequent substantial 
deposition of energy into the interstellar medium. 


To test this scenario we resimulated all our merger trees 


with the addition of BHs and with our implementation of 
BH feedback (see Sec. 12.211 . We also performed a simula¬ 
tion with BH but without the hot halo component for three 
merger trees, namely 780, 872, and 990. This is important 
to separate the effects of BHs and the hot halo and directly 
compare our results to previous studies. The top panels of 
Fig. [11] shows the results for merger tree 990 without the 
inclusion of the hot halo component. The upper left panel 
shows the SFR for three different simulations: host in isola¬ 
tion with BH, merger simulation with BH and merger simu¬ 
lation without BH. The SFR for tree 990 without a hot halo 
and without BH feedback is also shown for comparison (the 
same curve as shown in the left panel of Fig. 0. It is clear 
that the SFR rate is reduced in the presence of a BH, with 
the second peak lowered by more than a factor two. 

The effect of the BH is even more clear in the lower left 
panel, where the mass of the cold gas is shown as a function 
of time. With the inclusion of BHs and associated feedback, 
the cold mass decreases and quickly becomes negligible. The 
decrease of the SFR and the cold-gas mass can be easily 
linked to the increase of the gas-accretion rate onto the BH 
(top right panel in Fig. Hill . The accretion rate onto the BH 
is directly linked to the feedback energy output, responsible 
for heating the gas. T hese results are very similar to those of 
previous studies (e.g . Di Matteo et al.l l2005l : ISoringel et al.l 
I 2 OO 5 I : iHavward et al.ll2014l ~l. 

Next we included the presence of a hot gas halo in our 
initial conditions for merger 990 and reran it with BHs. The 
results are presented in the bottom row of Fig. 1111 The sim¬ 
ulation with BHs shows a similar trend as in the previous 
figure, when it is compared to the simulation without BHs 
but with hot halo. The total SFR is reduced and the second 
peak of SF (at 8 Gyr) drops by more than a factor 4. After 
the merger, SF is quenched to values as low as 2 Mq yr“^, 
while without a BH the SFR was around 10 Mq yr“^. 

The inclusion of the hot halo has a strong influence on 
the later stage of the evolution of SF and cold-gas mass. 
After 10 Gyr the hot halo is able to substantially increase 
the cold-gas mass in the center of the halo through cooling, 
as is deary shown in the top middle-left panel of Fig. 1111 
Thanks to this new cold gas reservoir, SF can increase again 
(contrary to the simulation with BH and no hot gas) and 
recovers to a SFR of several solar masses per year at the end 
of the simulation. 

The higher cold-gas fraction in the run with the hot halo 
does not only increase the SF burst after eight Gyrs (with 
respect to a run without the hot halo) but also increases 
the accretion rate onto the BH as shown in the lower right 
panel in Fig. 1111 This is true during the merger, but it is even 
more evident after (t > 9 Gyrs), when there is a difference 
of more than one order of magnitude in the accretion rate 
of the simulation with and without the hot gas component. 
These higher 7hi,h values result in a larger BH mass at the 
end of the run as shown in the bottom middle right panel 
of Fig. [TT] 

The effect of the BH on the gas distribution can be 
clearly seen in Fig. [T^l where the surface density map of the 
gas in the merger 990 is shown. Before the merger (t=6 and 
t=7 Gyr) the BH is depleting the gas content in the center of 
the galaxy, while during the merger (t=8 Gyr) a lot of gas is 
funnelled into the nucleus. This increases the accretion rate 
by about two orders of magnitude, releasing a lot of energy 
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Figure 9. Distribution of gas in tree 1166. The different snapshots represent the start of the simulation (t=4.64 Gyr), before, during, 
and after the multiple merger (t=7, 8 , 8.3, 9.3 Gyr), and the end of the simulation (t=13.7 Gyr). Different colours correspond to different 
gas surface densities as shown in the key. 
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Figure 10. The multiple mergers in trees 1188 and 1415. The leftmost columns show the SFR, the second column the total stellar 
mass, and the third column the cold-gas mass within a radius of 20 kpc. The rightmost column shows the average temperature T 30 . The 
red line is the central galaxy in the original simulation, and the green line shows the simulation where only the central galaxy and first 
satellite is included. The vertical orange dash-dotted and pink dashed lines (only in tree 1415) indicate the time when the two satellites 
merge with the host galaxy (/i = 1.36 for 1188, and p. = 0.87 and p = 0.69 in 1415) in the simulation with a triple merger, while the 
vertical dotted blue line corresponds to the time when the galaxies merged in the binary merger where only the most massive galaxies 
are included in the simulation (/i = 1.52 for 1188, and p. = 0.87 for 1415). 
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Figure 11. Merger tree 990 with BHs but without (top) and with a hot gaseous halo (bottom). For each row, the panels represent the 
following: in the left plots the SFR of the merger tree is shown, while the second plot from the left shows the mass of cold gas present 
in the galaxy. The right two panels show the BH mass and accretion rate. In the left two plots the solid red line shows the primary 
galaxy during the merging process, and the green line is the primary galaxy in isolation. Added to these plots are the results when BH 
feedback is not included, represented by the dotted red line. The red (blue) line in the right two plots shows the properties for the BH in 
the center of the primary (secondary) galaxy. When the galaxies and BHs merge the results for the merged BH continue as a single blue 
line. In the right plots we overplot a green dotted line, which corresponds to the ratio of SFR and mnh (merging galaxy), as indicated on 
the right vertical axis. The vertical lines in all panels indicate the time when the galaxies merge (orange dash-dotted for the simulation 
without BH, and vertical dashed purple for the simulation with a BH). 


from the BH that again heats and depletes the gas in the 
central region (t=13.7 Gyr). 

The reduced SFR in the simulations including BHs 
brings the final galaxies (i.e. the merger products) into bet¬ 
ter agreement with the observed stellar-halo mass relation 
as already shown in Fig.[^ in agreement with previous stud¬ 
ies that have shown AGN feedback is necessary to prevent 
the formati on of overly massive ga laxies, even at these mass 
scales (e.g. ISomerville et al.ll2008 ). 

Let us now summarise the results of the whole ensem¬ 
ble of simulated trees with BHs. The mass of the BHs in¬ 
creases slightly during the first stage of the simulation, the 
pre-merger phase. Lower mass BHs grow relatively faster 
than larger BHs as a consequence of using the Bondi Hoyle 
accretion parametrisation. 

The second growth phase is during the merger, when the 
accretion rate increases by an order of magnitude in most 
merger trees (see for example Fig. [13. We find a strong 
correlation between the presence and size of a starburst and 
the increase of the accretion rate of the BHs. Mergers that 
produce a large starburst also cause a large increase in the 
accretion rate, as both phenomena are related to the amount 
of cold gas in the galaxy center. 

Most of the BHs merge within the simulation time (this 
is visible for example in the middle right panels of Fig. 
m where one line disappears at around 8 Gyrs). After the 
merger the accretion rate is quite low even in the presence 
of a hot halo component and the mass of the BH remains 
roughly constant. 

These general trends are in good agreement with pre¬ 


vious studi es, e.g. Di Matteo e t al.l (l2005l ). ISoringel et al.l 
(l2005l) . and iHavward et al.l I 20l4) . with the difference that 
in their simulations SF is quenched indefinitely because they 
do not include a hot halo, while in our simulations we find 
that after several gigayears the halo has cooled enough to 
provide fuel for SF. 


The mass of the BHs at the end of the simulation can be 
compared to observations, for example using the correlation 
found between bulge mass and BH mass (e.g. Haring fc Rixl 
|2004 ISani et ^ I2OIII : iKormendv fc Hd l2013l ). In Fig. [TBI 
we show that there is reasonable agreement between the 
simulated galaxies and observations. This is due to the self- 
regulated n ature of the BH grow th as shown in previous 
studies (e.g. iDi Matteo et al.ll200^ . 


In two simulations with BHs, namely those of merger 
trees 1166 and 1684, we see that the above described be¬ 
haviour does not apply to every BH. In tree 1166, the third 
and smallest BH shows a fast and constant increase of its 
mass. In this merger tree there is an interaction between 
three galaxies and because of this, the third BH present in 
the smallest galaxy is able to move through the high-density 
gaseous medium of the other galaxies. This allows the BH to 
constantly accrete gas and keep increasing its mass, until it 
merges with the BH in the central galaxy. In tree 1684, the 
BH in the central galaxy shows barely any signatures of a 
merger and keeps increasing its mass exponentially through¬ 
out the entire simulation. This is probably due to the large 
amount of gas accreting onto the central galaxy. 


Finally, similar to the previous section, we want to 
quantify the ability of mergers in triggering additional SF 
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Figure 12. Gas surface density map of merger tree 990 with BH feedback included, shown at six different times. A depletion of gas 
can be seen around the BH (black empty circle) at all times, except during the merger at 8 Gyr, when more gas is transported into the 
nucleus of the galaxy, increasing the accretion rate of the BH. 


Tree nr 

Gyr 

t2 

Gyr 

h 

fgas 

f 

^bulge,host 

r 

•*-bulge,sat 

■^burst 

Mcold 
10^° Mq 

^burst 

661 

5.39 

6.32 

6.59 

0.63 

0.05 

0.61 

0.113 

1.875 

0.06 

661 

7.32 

10.72 

2.89 

0.38 

0.28 

0.00 

0.871 

3.491 

0.25 

872 

8.96 

9.35 

3.59 

0.20 

0.07 

0.62 

0.017 

1.820 

0.01 

990 

7.45 

8.50 

2.22 

0.17 

0.16 

0.24 

0.341 

1.175 

0.29 

1166 

7.18 

8.89 

3.30 

0.43 

0.05 

0.01 

0.184 

1.329 

0.14 

1166 

8.00 

8.89 

1.62 

0.27 

0.18 

0.03 

0.101 

0.854 

0.12 

1178 

7.40 

9.75 

2.59 

0.50 

0.02 

0.24 

0.666 

2.174 

0.31 

1188 

6.36 

8.11 

1.39 

0.42 

0.02 

0.04 

0.414 

1.466 

0.28 

1188 

10.27 

10.52 

8.14 

0.11 

0.85 

0.02 

0.002 

0.508 

0.004 

1188 

12.21 

13.72 

2.25 

0.12 

0.49 

0.17 

0.064 

0.604 

0.11 

1684 

8.45 

8.90 

2.32 

0.31 

0.11 

0.00 

0.028 

2.384 

0.01 

2096 

7.75 

10.70 

1.31 

0.36 

0.01 

0.002 

0.913 

2.284 

0.40 

3747 

5.20 

6.25 

4.35 

0.79 

0.01 

0.00 

0.285 

1.566 

0.18 

3747 

9.55 

11.50 

3.03 

0.23 

0.13 

0.11 

0.173 

2.011 

0.09 


Table 5. Simulations with BHs, only nine out of twelve merger trees showed major mergers. See caption TableO 


through the starburst parameter eburst- The relation be¬ 
tween the starburst efficiency and the merger ratio for sim¬ 
ulations with BHs is shown in Fig. 1141 while Table [5] is the 
analogue of Table [3] for these simulations. 

It is clear from the figure that the inclusion of both 
the hot halo component and BHs strongly reduces the effi¬ 
ciency of massive major mergers in t riggering additio nal SF 
with respect to previous studies (e.g ICox et al.ll2008l ). Both 
the slope 7 and the normalisation em of the power law are 
lower than found above, showing that the correlation be¬ 
tween starburst effiency and merger ratio is much weaker. 

The weaker correlation between the merger ratio and 
the burst intensity could be partially due to a self-regulation 
mechanism between SF and BH feedback. A massive merger 


is able to fun nel a significant amoun t of cold gas towards the 
nucleus (e.g. IPi Matteo et al ]H); this cold gas provides 
fuel for the stellar burst but is also available to be accreted 
by the BH. While in the absence of a central BH a very 
massive merger will trigger a strong star formation episode 
(as shown in figure [6|, in the presence of BH feedback the 
merger will also trigger accretion onto the black hole and 
the subsequent feedback, which will then suppress the stellar 
burst. 

Finally, the large scatter in the relationship between 
burst efficiency and merger mass ratio suggests that other 
factors can have an important role in determining the 
strength of the starburst. In Fig. [15] we again investigate 
if there is a single parameter that controls this scatter, but 
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Figure 13. Bulge mass versus BH mass for the simulated galaxies 
at z = 0 (blue points). The black lin e shows the observed corre¬ 
lation found by iHaring &: Rii3 ll2004h . while the grey area shows 
the intrinsic scatter in this relation. 



Figure 14. Merger mass ratio versus starburst efficiency. The 
blue points show all major merger events in simulations with a 
BH and a hot gaseous halo around the galaxy, while the green 
squares show the simulations with black holes and withou t a hot 
halo. The black line is the starburst efficiency of ICox et al.l (l2008l ) 
and the green line represents the least squares fit to our results. 


Tree 

tl 

Gyr 

t2 

Gyr 


■l^burst 

Mcold 
10 i° Mq 

^burst 

872 

10.95 

12.00 

3.39 

0.001 

0.298 

0.000 

990 

7.30 

9.30 

2.38 

0.220 

0.699 

0.315 


Table 6. Results of the simulations without a hot halo, but with 
a black hole. Only two out of twelve merger trees were simulated 
with these initial conditions. The columns correspond to columns 
1, 2, 3, 4, 8, 9, and 10 from table|3] 


we only find a possible weak trend with At, an indicator of 
the orbit. We plan to address the origin of this scatter in 
more detail in a future paper. 


6 DISCUSSION 

This study differs from previous ones in our method of ob¬ 
taining initial conditions (IC) for numerical hydrodynamic 
merger simul ations. Previous studies that included a hot 
gaseo us halo JCox et al.ll2006l : lMoster et al.l[20lI] : IChoi et al.l 
I2OI4I ) had ad hoc ICs, while here we adopted ICs that are 
cosmologically motivated. Therefore, the progenitors of our 
galaxies have a broad range of initial properties (gas frac¬ 
tion, morphology, size, etc.), and consequently, our results 
are more representa t ive of an unbiased galaxy population. 

iHavward et al.l ll2014l) compared the performance of 
SPH to AREPO in merger simulations. They showed that 
SPH poorly resolves thermal instabilities, and it produces 
a signifi cantly larger h ot gaseous halo around galaxies than 
AREPO l|Springelll20ld ). with spurious blobs of cold gas. The 
gravitational heating by satellites that we see in our simu¬ 
lations might be sensitive to the these numerical issues, and 
should be more thouroughly tested with grid-based codes 
or improved versions of SPH. In addition, the blobs formed 
in SPH might spuriously contribute to the heating of the 



bulge 


At 


Figure 15. Difference between the burst efficiency compared to 
the power law fit to the data in Fig. 1141 Aeturst = Cfaurst ~ 
Cburst.PL- The difference is shown as a function of the gas fraction 
/gas (left), the bulge fraction of the main galaxy /bulge (middle), 
and the duration of the starburst At = t 2 — ti which is a tracer 
of the orbit (right). The big circled points are all major mergers 
with 1 < /r < 3, while the smaller points are all mergers outside 
of this range. 


hot halo, although it is likely a minor effect, as most of the 
heating occurs shortly before a merging event. 

Although we find clear evidence that the recent merger 
history should be taken into account when calculating burst 
efficiencies, a full study of the effects of multiple mergers 
would require an additional substantial simulation effort, 
which is beyond the scope of this paper. We plan to address 
this in a forthcoming work. 

In our suite of simulations, we find some significant star- 
bursts for relatively large merger ratios (p > 5, see Fig. HI, 
contrary to what previous studies found. Because we focused 
on major mergers in this study, our selection criteria limit 
the number of minor mergers. Consequently, we do not have 
a large enough sample of minor mergers to constrain the 
dispersion in burst efficiency for minor mergers. 

In Figs. [5]and|8l we see that the gas mass for the galax¬ 
ies is nearly constant with time. Because the stellar mass 
does increase by a factor of ~ 2 since z = 1, the gas frac¬ 
tions of these galaxies also decrease by a factor of 2. Since we 
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are looking at individual galaxies, it is difficult to compare 
the gas masses to observations. However, the trend that the 
gas fraction decreases by a factor of 2 over this per iod is 
compatible with observations fe.g. ISantini et al]l2014h . The 
constant cold gas masses are also in agreement with observa¬ 
tions that the neutral-gas content of galaxies does not vary 
in ga laxies that recently underwent a merger (lEllison et al.l 

l20ldl . 

From the stellar mass to halo mass comparison, we find 
that overcooling is not a major problem for the merger tree 
simulations discussed in this paper. However, this is in part 
due to the setup of our simulations, because they are based 
on SAMs that do include radio mode feedback. In addition, 
the rather large value of a, the ratio of specific angular mo¬ 
mentum of the hot halo gas and the dark matter, was cali¬ 
brated to obtain galaxies with reasonable star formation his¬ 
tories fsee lMoster et aHl201lh . For a test simulation with a 
significantly larger halo mass, we find very high SFRs caused 
by overcooling, indicating that at higher masses additional 
feedback such as radio mode AGN feedback is required (e.g. 
ISomerville et al.]l2008l : iFanidakis et al.ll201ll l. 

We found, in agreement with previous studies, that the 
inclusion of AGN feedback reduces the efficiency of merger- 
triggered bursts. In our simulations, black hole accretion and 
the associated feedback are modelled in a rel atively crude 
manner. Recent studies (e.g. IChoi et al.ll2014h have shown 
that the deposition of momentum from AGN-driven winds 
originating on the unresolved scales of the nucleus may have 
an even more dramatic effect, driving gas out of the nu¬ 
cleus in high-velocity outflows and further suppressing nu¬ 
clear star formation. 


7 CONCLUSIONS 

Starbursts triggered by galaxy mergers have long been a 
topic of great interest in the literature. Here, we present 
important progress in the theoretical understanding of this 
phenomenon by improving on past studies in several ways. 
In particular, we employed the Simulated Merger Tree ap¬ 
proach introduced in MMS14, in which cosmological merger 
trees combined with a semi-analytic model are used to ob¬ 
tain the initial conditions (orbit and galaxy properties) for 
a sequence of hydrodynamic merger simulations carried out 
with the SPH code GADGEt2. As motivated by fully cosmo¬ 
logical simulations and observations, but unlike most pre¬ 
vious studies, our model galaxy-halo systems include a hot 
gaseous halo. We studied the impact of the inclusion of the 
hot gas halo, as well as that of BH feedback and multiple 
mergers, on triggered star formation enhancements in galaxy 
interactions, with a focus on major mergers. 

Our main results can be summarised as follows: 

• Owing to the presence of the physically and observa- 
tionally motivated hot halo, our simulated galaxies are able 
to sustain SF for longer times, without experiencing artifi¬ 
cial strangulation due to lack of cold gas. 

• The starburst which results from the merger of two 
galaxies is relatively weak compared with most previous sim¬ 
ulations. This is mostly because of the increase of the SFR 
at later times (due to the hot halo), which does not decrease 
the absolute size of the starburst, but does decrease its rel¬ 
ative size. 


• In some of our mergers the total SF is decreased with re¬ 
spect to the SF of the same galaxy run in isolation. This im¬ 
plies that mergers can also have a negative effect on SF. This 
decrease of SF is caused by the transfer of orbital energy of 
the merging galaxies into thermal and kinetic energy of the 
hot gas and by the generation of shocks in the hot halo dur¬ 
ing the merger. These heating processes increase the cooling 
time of the hot circumgalactic gas and lead to a lower SFR 
in the latest stage of the mergers. This also shows that ex¬ 
tra heating terms are present during merging processes and 
should be taken into account in galaxy evolution models. 

• We showed that three simultaneously merging galaxies 
should not be treated as two independent and consecutive 
mergers. The starburst resulting from a triple merger can be 
either similar to the one of the merger of the two most mas¬ 
sive galaxies or can even be smaller. These results underline 
once more the complex dynamics of multiple mergers and 
the non-linear effect of different heating and cooling mech¬ 
anisms. 

• When we include thermal feedback from massive cen¬ 
tral BHs, we find that the efficiency of merger-triggered star- 
bursts is reduced, especially right after the merger. Overall, 
BHs substantially reduced the stellar mass of the merger 
remnants, bringing them in better agreement with abun¬ 
dance matching predictions. 

• Finally, following previous studies, we described the ef¬ 
ficiency of mergers in triggering a starburst by a power-law 
function of the merger mass ratio. In absence of BHs, we 
find good agreement in the normalisation of this relation 
with iMoster et ^ ll201lll who found em = 0.51 in their 
hot halo simulation, while we find em = 0.52 for a larger 
sample of galaxies. Our value for the index of the power-law 
7 = 0.81 is in rea sonable agreement with the value found by 
ICox et ^ ll2008ll (7 = 0.69). However, the scatter that we 
find is significantly increased, arguing against using a sin¬ 
gle fit to determine burst efficiencies for individual galaxies. 
When we include feedback from BHs there seems to be a 
weaker dependence of the starburst efficiency on the merger 
mass ratio. 

There are still numerous improvements possible for our 
simulations. For example, we have a quite crude description 
of the thermal feedback from the central black hole and we 
did no t include mome ntum deposition f rom nuclear scale 
winds llChoi et al.ll2014ll or any jet-like fe.g ICielo et al.l[20l3 l 
or radio feedback from the BH. 

Nevertheless our results clearly show the importance of 
including all galaxy components (DM, stars, hot gas, cold 
gas, BH) when performing galaxy merger simulations, and 
of adopting cosmologically motivated merger histories and 
initial conditions. Our results may be useful for improving 
the treatment of merger-driven starbursts in semi-analytic 
models of galaxy formation. 
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APPENDIX A: NUMBER OF PARTICLES IN 
SIMULATIONS 
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Tree nr 

type 

9 ^ gas particles 

# DM particles 

9 ^ stellar disc particles 

9 ^ stellar bulge particles 
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h 
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334836 
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0 
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Table Al. The number of particles used to generate galaxies in our merger trees. Central galaxies are labelled in the second column 
with h, while the satellites are represented by s. 



